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^ ■ Abstract 

. A recent analytic test of the instanton method performed by comparing the 

O ! exact spectrum of the Lame potential (derived from representations of a 

I finite dimensional matrix expressed in terms of su{2) generators) with the 

O ' results of the tight-binding and instanton approximations as well as the 

^ standard WKB approximation is commented upon. It is pointed out that in 

. the case of the Lame potential as well as others the WKB-related method of 

! matched asymptotic expansions yields the exact instanton result as a result 

OO ! of boundary conditions imposed on wave functions which are matched in 



o 



domains of overlap. 



Oh! 
(D . 

-F; ■ 1 Introduction 

^ ! We comment on ref.|l|. Ever since the instanton method of evaluating path 

H ' integrals became popular, the question was raised as to whether this method 

and that of the Schrodinger equation lead to identical results in the domi- 
nant approximation. In view of the significance of both methods this is an 
important question. Generally one would argue that in the 1-loop approxi- 
mation of the path integral the results of both methods should agree - and, 
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in fact, this is the widely accepted opinion. However, there are few mod- 
els which allow explicit verification, and some detailed investigations of even 
these are not well known. Thus recently in a study of Lame instantons [Q], 
a parallel consideration of the standard WKB method (i.e. that with linear 
connection relations) was found to lead to a result which is off by an overall 
factor of ^e/vr - very similar to an observation made long ago in ref . [0] in the 
case of the well known double well potential. In the following we point out 
that in both cases, the WKB method employed is too simple, and complete 
agreement between the result of the path integral method and that of the 
Schrodinger equation in all the examples referred to can be achieved with 
the method of matched asymptotic expansions Q 

The best known examples for such considerations are the scalar theory 
with the double well potential and the sine-Gordon theory of the periodic 
potential. The case considered in ref.|I| with the Lame potential is not 
so well known but has recently found widespread consideration in various 
contexts, such as supersymmetric quantum mechanics |^ and spin tunneling 
(cf. e.g. P), and includes in a particular limit the sine-Gordon case. In 
the following we sketch the main points of the method and recall the level 
splitting calculated long ago for the Lame (and more generally ellipsoidal) 
wave equation |0] even for an arbitrary excited state, and consider limits and 
special cases to demonstrate the agreement with the result of the instanton 
method, or more generally the path integral method, in particular that of 
the ground state case of ref.[|l|. For details of the calculations we refer to the 
literature, where many details are given. 



2 The Lame level splitting and consequences 

The Lame equation p 

y" + |a - K^sn^u^y = 0, = n{n + 1)A;^ (1) 

*The source of the discrepancy can be traced to a mismatch between the normahzation 
of the WKB wave functions and the harmonic oscillator wave functions, as explained in 
ref. and also commented on in ref. We thank the referee for pointing this out. 
Our solutions involving the function i3[z(ti)] in terms of Hermite functions correspond to 
the harmonic oscillator wave functions of ref. . Matching WKB- or WKB-like functions 
(like our functions involving A{u)) to these therefore avoids the mismatch 
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with elliptic modulus \k\ < 1 and n > —1/2,0 < u < 2/C, can be looked 
at as a Sclirodinger equation with periodic potential n'^sn'^u, where snu is 
one of the Jacobian elliptic functions with period 2/C, and /C is one of the 
complete elliptic integrals. (In the comparison with the Schrodinger equation, 
the usual factor —h'^ /2m in front of the second derivative has to be kept in 
mind, where m is the mass). 

We recapitulate briefly the main steps in the method of matched asymp- 
totic expansions (of wave functions) as applied to the periodic Lame equation 
(IID (for the simplest application of the method see ref. 0). The first step is 
to write the eigenvalue A as 

A(g,K):=g«: + ^^, (2) 

where q ^ Qo = 2N + 1, = 0, 1, 2, ■ ■ ■ in the case of k ^ oo, i.e. very high 
potential barriers (i.e. harmonic oscillator approximation around a minimum 
of the potential). For barriers of finite height the parameter q is only approx- 
imately an odd integer go in view of tunneling effects. The difference q — qo is 
obtained by imposing boundary conditions at extrema of the potential (see 
below). The second step is to insert @) into the equation, i.e. (P, and to 
write 

y = y4(u)exp| — J Ksnudu^ 

= A{u)[fiu)r/''^ (3) 

where 

^ /dnu + kcnu\ 
\dnu — kcnuj 

For large k the equation for A{u) can be solved iteratively resulting in an 
asymptotic expansion for A{u) and concurrently one for the remainder in eq. 
(H), i.e. A. A second solution is written 

y = A(u)e.p{+lnsnudu}. (4) 

The very useful property of these solutions is that for the same value of A 
(which remains unchanged under the combined replacements q —q, k, 
-k) 

A{u) = A{u + 2}C), 
A{u,q,K) = A{u,-q,-K). (5) 
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The domain of validity of these solutions is that away from an extremum of 
the potential, more precisely for 



dnu =p cnu 



dnu ± cnu 



1 

>> -. 

K 



Thus one can construct solutions Ec{u), Es{u), which are respectively even 
in u (or snu) or odd, i.e. 

^A{u)[f{u)r/''±A{u)[f{u)]-^/'' (6) 

Since these expansions are not valid at the extrema of the potential (where 
the boundary conditions are to be imposed), one has to derive new sets of 
solutions there and match these to the former (i.e. determine their propor- 
tionality factors) in domains of overlap (their extreme regions of validity). 

Thus in the third step two more pairs of solutions B, B and C, C replacing 
A, A are derived, one pair in terms of Hermite functions of a real variable, 
the other in terms of those of an imaginary variable, by transforming the 
equations for A, A into equations in terms of 

k \ dnu ± cnu J 

Solving the resulting equations itcratively as before, one obtains again the 
same expansion for A, but solutions B, C replacing A, which are valid for 



dnu ± cnu 



dnu ^ cnu 



« 1. 



In their regions of overlap one can determine the proportionality factors a, a 
of 

B = aA, C = 6lA (8) 
(again in the form of asymptotic expansions). Then 
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In the fourth and final step one now apphes the appropriate boundary 
conditions (cf. ref. p[) on these solutions (j^), i.e. one sets 

Ec{u = 2/C) = Ec{u = 0) = 0, 
Es{u = 2/C) = Es{u = 0) = 0, 



dEc\ _ /dEc\ 

du ) u=2K. \ du ) u= 



0, 



( du )«=2c ( du )„=o ^' ^^^^ 

These conditions define respectively functions of period 4/C, 2/C, 2/C and 4/C. 
Evaluating these one obtains (from factors of factorials in q and —q) expres- 
sions cot{7r(g — l)/4} = ■ ■ ■ , tan{7r(g — 1)/4} = ■ ■ ■, from which the difference 
g — go is obtained by expansion around zeros. Finally expanding 



A(g) ~ A(go) + (g - go) 



{dq) 
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one obtains the eigenvalues from which the level splitting can be deduced. 

In view of its periodicity the periodic potential avoids the necessity of 
matching across turning points in the above calculation. In the case of the 
double well potential this is different, as explained in detail in ref. [|T^], 
and one has to impose boundary conditions not only at the minimum but 
also at the central maximum, again, of course, on even and odd solutions 
constructed parallel to those above. 

The perturbatively derived wave functions of ref.0 (for large values of 
K^), when matched in domains of overlap and so extended over the entire 
domain of the variable m, and subjected to periodic boundary conditions as 
described above define two pairs of eigenfunctions, in each case with one 
even and one odd, of periods 2/C and 4/C respectively. These four conditions 
together imply for large values of the following asymptotic expansion of 
the eigenvalues A as shown in ref. [0: 



A±(go) =A(go)±TY- 



l + k\-kf 8k \29o 



kj \l-k 



l + Ol- 



K 



Here go = 2A^-|- 1, = 0, 1, 2, ■ ■ ■ and A(go) is the purely perturbative contri- 
bution which represents effectively the eigenvalues of degenerate oscillators 



5 



of the periodic potential in the case of very high barriers. It is the boundary 
conditions imposed on the perturbatively derived solutions which yield the 
nonperturbative effects equivalent to those of the instanton. Thus the factor 

1 + k\-'^/^ 



1 - k 

is, in fact ea;p(— S'o), where 5*0 is the Euclidean action of the instanton. 

We first verify the result (|TT]) by reduction to the sine-Gordon case. With 
K = ±2h finite while n ^ oo and k ^ the Jacobian elliptic function snu 
reduces to sinw and eq.(|l]) becomes by replacing u hj x ± n/2 the Mathieu 
equation 

y" + S^\-2h^cos^2x^y = 0, \ = A-2h^. (12) 
Under the conditions stated the eigenvalues become 



A±(go) = A(go) ± 



(13) 



in agreement with established results in this case |1^, p!l[| . 

Next we consider k approaching 1 in the case of the Lame eigenvalues 
([lID (terms up to and including those of 0(1/k^) in the level splitting and up 
to and including those of 0(1/k^) in the perturbative part have been given 
in ref. for any go)- One readily obtains 



A±(go) = A(go)±7Y- 



5/t 2 



+1/ 



go 



L2-(go-l)]!(27r)i/22«+i+f 
For the two lowest levels qo = I and one obtains 



l + (l-A;){«:(l-ln2) + |}+0 



I 

(14) 



A±(1) = A(1)± 



(4/€)3/2(l - kY~-2 



(27r)i/22 

Thus the separation of the two lowest levels is 
2(4k)3/2(1 - k)^-^ 



l + (l-fc){fi:(--ln2) + -}+0 - 



(15) 



AA(1) 



(27r: 



l + (l-A:)wl-ln2) + l} + 0(i) 



(16) 
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This result agrees with formula (13) of ref.|I| in the limit of (in our notation) 
large k and /c — 1 (in particular this agrees also with our power k— l/2of(l — 
fc) where in the —1/2 has been ignored). The higher order contributions 
are, presumably, somewhat model or approximation dependent. We thus 
have agreement with the instanton result of ref. fl]] (there v is our fc^, so that 
for /c — > 1 one has (1 — z/) — > 2(1 — fc)). 

Hence in the case of the periodic cosine potential, as well as in the case of 
the Lame potential, the method of matched asymptotic solutions of refs. [0, 
- yields the same result as the instanton method in the 1-loop approximation, 
as one would expect. In ref.[|I| reference is made to the well known case of the 



double well potential. For this case also it has been shown in ref. [[T^] that 
the method of matched perturbation expansions - which one might argue 
amounts to an improved form of the standard WKB method (since it gives 
the correct eigenvalues) - yields the same result as the instanton method. 



again for any arbitrary level, whether ground state or excited |]I0 



3 Conclusions 

In the above we have demonstrated that the method of matched asymptotic 
expansions of refs. 0, ||, |lO| which, incidentally, was also developed and used 
to determine the large order behaviour of the perturbation expansion (cf. 
ref.||l3|), leads to the same result as the instanton calculation in the 1-loop 
approximation and thus is superior to a WKB calculation. The essential 
difference is that WKB solutions alone, or rather our WKB-like solutions in- 
volving the function A, (which are not valid at an extremum of the potential) 
do not suffice; one has to match these to solutions valid around the extrema. 
We may conclude that to obtain from the Schrodinger equation and pertur- 
bation theory results agreeing with those of the path integral method with 
expansion around the instanton, one has to use the full method of matched 
asymptotic expansions as developed in refs.0 and [|^ |T0[ (which has also 
been applied to other cases such as spheroidal wave equations [|l4l). In purely 
quantum mechanical cases, such as those considered above, the method of 
matched asymptotic expansions seems to be simpler and yields the splitting 
of excited oscillator states with the same ease as that of the ground state, 
whereas in the case of the path integral method, one has to use periodic 
instantons, as in e.g. refs. |1^, |T5|. This may be worth noting in connection 
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with models of spin tunneling which attracted considerable interest recently, 
since these can - in certain cases and with certain approximations (and co- 
herent states) - be related to periodic differential equations. Thus in ref. 
the case of the Hamiltonian H = KiS^ — K2SI describing a ferromagnetic 
particle with large spin has been considered and related to Mathieu and 
Lame equations. In particular the specific Lame instanton of ref. [|I| has 
been used in ref.|]l6|. Finally we add that the introduction of the parameter 
q in conjunction with the construction of solutions with the properties of eq. 
@ has been shown to be extremely useful in other but related contexts as 
may be seen from the calculation of the scattering matrix in ref.|]l7 . 
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